Anomalous Finite Size Effects on Surface States in the Topological Insulator Bi 2 Se 3 
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We study how the surface states in the strong topological insulator Bi2Se 3 are influenced by finite size effects, 
and compare our results with those recently obtained for 2D topological insulator HgTe. We demonstrate two 
important distinctions: (i) contrary to HgTe, the surface-states in Bi2Se3 display a remarkable robustness to- 
wards decreasing the width L down to a few nm, thus ensuring that the topological surface states remain intact, 
and ( ii) the gapping due to the hybridization of the surface states features an oscillating exponential decay as a 
function of L in Bi2Se3 in sharp contrast to HgTe. Our findings suggest that Bi2Se3 is suitable for nanoscale 
applications in quantum computing or spintronics. Also, we propose a way to experimentally detect both of the 
predicted effects. 

PACS numbers: 73.43.-f, 72.25.Dc, 85.75.-d 



I. INTRODUCTION 

A new state of matter, known as the topologically insulating 
state, has recently been experimentally observed^ 3 -! after the 
successful prediction of its existence in HgT#l This state is 
characterized by the topological protection of the conducting 
states that form at the edges (in 2Dj^C or surfaces (in 3Dj2U 
of such materials, whereas the bulk states remain insulating 
due to a charge excitation gap. This distinguishes topologi- 
cal insulators from conventional insulators that do not feature 
such e dge/surfa ce states, and introduces several interesting 
effects 10 ! 11 ! 12 ! 13 ! . In particular, topological insulators could find 
use in quantum computation since the topological ly protec ted 
edge/surface states remain insensitive to disorde r 4 ! 15 ! 16 ! 17 !. 

It has recently been realized that Bi 2 Se 3 is a three- 
dimensional topological insulator with a large charge exci- 
tation gap in the bulk 18 . The surface states have an energy 
dispersion that is linear in momentum and thus form a Dirac 
cone at low energy, similarly to graphene. In stark contrast, 
however, Kramer's theorem does not guarantee the survival of 
the edge states in graphene in the presence of perturbations 
since it holds an even number of Dirac cones inside the Fermi 
contour. 8 In Bi2Se3, the number of Dirac points inside the 
Fermi arc is odd which activates the protection from Kramer's 
degeneracy theorem. This has very recently been experimen- 
tally observed-51. 

The prospect of utilizing the protected surface states in 
topological insulators such as Bi2Se3 in actual devices re- 
lated to quantum computing or spintronics demands that finite 
size effects are taken seriously. This fact is underlined by the 
finding of RefPS which showed that the edge states of HgTe 
quantum wells become gapped due to a finite-size effect as 
the width L is decreased. The gapping was shown to become 
experimentally measurable around L ~ 200 nm, suggesting 
that the material loses its exotic edge state properties in this 
region. Clearly, this places severe restrictions on potential use 
of HgTe quantum wells in applications on the nm-scale. 

In this work, we will demonstrate how the situation changes 
dramatically when considering Bi2Se3. Within a combined 
analytical and numerical approach, we show how the sur- 
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FIG. 1: (Color online) Suggested experimental setup: two-terminal 
geometry with the topological insulator Bi2Se3. The direction of the 
current carried by the surface states is shown by the arrows, and the 
width of the sample is L. For L ~ 2hvF /\M\, finite-size effects 
become important. Here, \M\ is the bulk charge excitation gap. 

face states in Bi2Se3 display a remarkable robustness towards 
finite-size effects, becoming measurably gapped only when 
the width drops to a few nm. This means that samples of 
Bi2Se3 can be made several tens of times smaller than HgTe 
while still retaining their characteristic surface states giving 
rise to the quantum spin Hall effect. We explain this observa- 
tion in terms of the large charge excitation gap \M\ in Bi 2 Se 3 
which gives rise to a very short localization length of the sur- 
face states. Moreover, we show how the finite-size induced 
gapping in the surface states displays a qualitatively differ- 
ent dependence on L in Bi 2 Se 3 compared to HgTe, namely 
an oscillating exponential decay. Both of these effects can be 
measured in a two-terminal geometry as sketched in Fig. [T] 

II. THEORY 

The effective low-energy Hamiltonian for Bi2Se3 centered 
around the T point in the Brillouin Zone may be written a^l 

^ _ f£kl + M k Tz + AxkzTx A 2 k-Tx \ 

(1) 
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FIG. 2: (Color online) Plot of the surface/edge state energy disper- 
sions versus the transverse momentum k± in the case of Bi2Se3 and 
HgTe. For Bi2Se3, we give only the result for L = 1000 nm since 
the dispersion remains completely unchanged down to L — 150 nm, 
in contrast to HgTe where a gap opens up. The inset of Bi2Se3 shows 
the bulk bands and the surface states with a gap of 2\M\ = 0.56 eV 
opening at the V point. 



where we have defined the following quantities: 



e k = C + D 1 k 2 z + D 2 k 2 ± ,, 
M h = M - B x k 2 z - B 2 k\, 



(2) 



and k± = k x ± ik y . Here, r 2 and t x are the Pauli matrices in 
standard notation, while the parameters {Aj , Bj , C, Dj , M} 
describe the band-structure in Bi2Se3, obtainable by first- 
principles calculations. Such a fitting procedure was un- 
dertaken in Ref. 18 , with the result C = -6.8 x 10~ 3 eV, 



M = 0.28 eV, and 

A 1 = 2.2 eVA, A 2 = 4.1 eV-A, B l = 10 eV-A 2 , 

B 2 = 56.6 eVA 2 , D x = 1.3 eVA 2 , D 2 = 19.6 eV-A 2 . 



(3) 



It should be noted that there is an anisotropy along the i:-axis 
and that the full 3D structure of Bi2Se3 has been taken into 
account. The basis ^ we have used for the 4x4 Hamiltonian 
is 

* = (|P1+ t), \P2~, T), |P1+ |), \P2~, j)) tr , (4) 

where tr denotes the transpose operation. The above p-orbital 
states are the relevant ones near the Fermi level of Bi2Se3 and 
may be classified according to their parity ± since inversion 
symmetry is preserved on the lattice. The states |P1) stem 
from the Bi atoms, whereas \P2) stems from the Se atoms. 
The |P1 + ,(t) and \P2~,<r) states (a =t,4) have opposite 
parity and the order of them near the Fermi level is inter- 
changed when spin-orbit coupling is taken into account. In 
this way, the spin-orbit coupling effect is responsible for driv- 
ing the system into a topologically insulating phase. 

It is interesting to observe that the Hamiltonian Eq. ([T]) 
contains spin-mixing terms due to the off-diagonal entries 
A 2 k±r x . This in contrast to the 2D topological insulator 
HgTe, which is diagonal in spin-space. As a result, one might 
expect new features in the spin-current of Bi2Se3 carried by 
the topological surface states. 



III. RESULTS AND DISCUSSION 

We will consider a finite width L in the z-direction with 
open boundary conditions at the edges, i.e. 



%(x,y,z= ±L/2) = 0. 



(5) 



Since the translational symmetry is broken along the z- 
direction, we perform a Peierls substitution k z — * — \d z in 
the Schrodinger equation = e^. Assuming a plane-wave 
solution 



(6) 



where A determines whether the mode is evanescent or prop- 
agating, we solve the secular equation to find the allowed 
eigenvalues for A. These read A = ±A Q , where 



A Q = [A\ - 2£>i(C - e + D 2 k\) + 2B x (B 2 k\ - M) 



+ a^R] l l 2 /^2(Bl-Dl), 
with the definition 



±1, 



(7) 



R = Aj[Af -AD 1 {C~e + D 2 k\) - AB^M - B 2 k 2 ± )} 
- AA\k\(B\ - D\) + A[B X {C - e + D 2 k\) 
+ D 1 {M-B 2 k 2 A _)] 2 . (8) 
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As demanded by consistency, Eq. (|7]i reduces to the result 
pl j n me limiting case of zero anisotropy and C = 0. The 
total wavefunction ^ is a superposition of the terms e ±AaZ 
with belonging normalization coefficients. The open bound- 
ary conditions at z = ±L/2 allow us to write down an im- 
plicit equation for the energy eigenvalues of ^. The values 
of e solving this equation then correspond to the bulk states 
and, where possible, surface-bound states. We arrive at the 
following energy eigenvalue equation: 

^ ta _ h\+^-(k\-^Y(Bl-Dl)/Al 

a 

where t a = tanh(A a L/2). For arbitrary values of k±, Eq. (|9ji 
cannot be solved analytically. Instead, we employ a numerical 
solution of Eq. (|9]l to find the allowed values of e for a given 
value of k±. All other material parameters are specified in Eq. 

To highlight the dramatic difference between the surface 
states in Bi2Se3 and the edge states in HgTe, we show in Fig. 
[2] the energy dispersion as a function of the transverse mo- 
mentum k±. The plots for HgTe were obtained by utilizing 
the results in Ref. 20 . In the lower panel, it is seen how a gap 
A opens between the edge states at the V point k± — as L 
decreases. Let us emphasize here that the gap A between the 
edge/surface states is to be distinguished from the gap \M\ 
between the bulk energy bands. For large L ~ 1000 nm, 
the edge states remain ungapped for HgTe. In sharp contrast, 
the edge states in Bi2Se3 remain completely ungapped in the 
entire regime L £ [150, 1000] nm considered in Fig. [5] In 
fact, we find that a measurable gap A does not begin to open 
until widths of L < 10 nm are reached. This fact suggests 
that much smaller samples that retain their conducting surface 
states can be fabricated in the case of Bi 2 Se 3 than in the case 
of HgTe, which is our first main result. Such an observation 
is crucial for the prospect of utilizing topological insulators 
such as Bi2Se3 and HgTe in applications linked to quantum 
computing or spintronics. 

In order to understand the large quantitative difference be- 
tween the necessary width L that induces gapping between 
the surface/edge states in Bi2Se3 and HgTe, we study more 
carefully the eigenvalues A Q . The physical interpretation of 
these quantities is that Re{A Q } corresponds to an inverse lo- 
calization length (or, alternatively, penetration depth into the 
bulk) for the surface states. Therefore, the largest of the length 
scales (Re{A + }) _1 and (Re{A_}) _1 mainly determines the 
density profile for the surface states and their penetratation 
into the bulk. We have verified numerically that 

(Re^})- 1 > (RejA+l)- 1 (10) 

for all energies inside the bulk gap ±|M| for both Bi2Se3 and 
HgTe, so that a = — 1 will determine the penetration depth of 
the surface states into the bulk. The penetration depth of the 
surface states can be estimated by 

£ = hv F /\M\, (11) 

where the Fermi velocity is provided by Vp — A/h for HgTe 
and v F = A 2 /h for Bi 2 Se 3 . 218 . Since A ~ A 2 while 
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FIG. 3: (Color online) Plot of the finite-size induced gap A between 
the surface states at k± = for Bi2Se3 as a function of the width L. 
In contrast to HgTe where the decay with L is purely exponential, an 
oscillatory pattern is superimposed on the decay for Bi2Se3. 



|-W|Bi 2 Se 3 3* |-W|HgTe> the large difference in the distribu- 
tion length of the surface/edge states stems from the sizable 
charge excitation gap \M\ in Bi2Se3. It is tempting to draw an 
analogy to the midgap Andreev-bound states in <i-wave super- 
conductors induced at the interface which extend a distance 
into the bulk proportional to the inverse of the superconduct- 
ing garP. 

Our second main result is related to the manner in which 
the gap between the surface states in Bi2Se3 depends on the 
width L, in effect A = A(L). It is instructive to first recall 
that A(L) in HgTe was shown to exhibit a purely exponen- 
tial decay in Ref. 20 . To investigate the situation in Bi 2 Se3, 
we have solved numerically for A(L) and plotted the result in 
Fig. [3] It is seen that a qualitatively different scenario from 
HgTe transpires: the decay with L is highly non-monotonous, 
and in fact features a superimposed oscillatory pattern on the 
exponential decay. The experimental signature of such an os- 
cillatory decay would be to measure the conductance at a fixed 
chemical potential for several samples with different widths L 
and see how the conductance appears and reappears, as we 
shall describe below. 

We now proceed to explain the origin of the oscillatory 
decay of the gap A(L) found in Bi 2 Se3, considering the T 
point k± = 0. The crucial observation in this context is that 
the eigenvalues A Q are not purely real in the bulk insulating 
regime e S [C — \M\, C + \M\]. This is in contrast to HgTe, 
where A Q are purely real in this regime. For Bi2Se3, we find 
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FIG. 4: (Color online) Plot of the two-terminal conductance at T = 20 K and T = 200 K as a function of /i and L. 



that 



A + = Alforee [C- |Af|,<7+ |M|], 



(12) 



as can be verified directly from Eqs. ([3]l and (|7]i since i? < 0. 
As a consequence, whereas a gap dependence of the type 
e found for Hg dictates an exponential decay, it will 
give a superimposed oscillatory pattern on top of the exponen- 
tial decay in the case of Bi 2 Se 3 . The natural question is then: 
is it possible to identify the reason for why R < in Bi2Se3 
whereas R > in HgTe, leading to qualitatively different be- 
havior of the surface states gap? Analyzing the expression for 
R in Eq. ([7]> at /cj^ = 0, it seen that neither the 3D nature or 
the anisotropy of the former material can be the reason, since 
there is no mixing between indices with subscript '1' and '2'. 
Therefore, R < seems to occur as a direct result of the ma- 
terial parameters given in Eq. rf3j. In principle, since all the 
parameters in the effective Hamiltonian of HgTe depend on 
the thickness d of the sample (which also determines whether 
the material is in the trivial or topologically insulating state), 
it might be possible to obtain a conversion from exponential 
decay to oscillating decay also in HgTe, although this clearly 
warrants a separate investigation based on first-principles cal- 
culations. 

The simplest way to experimentally detect both the gapping 
of the surface states and their unusual dependence on L is 
arguably a two-terminal geometry, as shown in Fig. [T] In 
that case, the conductance can be evaluated in the Landauer- 
Biittiker framework at a finite temperature T. Considering a 
zero-bias situation eV — ► 0, the transmission coefficient T 
may be written as 

T(e) = N c [e(e - A/2) + 6(-e - A/2)], (13) 



where is the Heaviside-step function, N c is the number of 
conducting channels on the surfaces z — ±L/2, and e is the 
quasiparticle energy. Similarly to RefP, one arrives at the 
following expression for the conductance normalized against 



G = -{I + [1 + c /3(A/2-M)]-l _ [! + ^(-A/2-p)]-l }) 



(14) 



with /3 _1 = ksT. As a direct consequence of the gap A, the 
conductance is suppressed as T — > 0. In Fig. |4] we plot the 
conductance as a function of L and the chemical potential /i, 
comparing the temperatures T = 20 K and T = 200 K. As 
seen, the conductance displays oscillations for a fixed \i, and 
the oscillation length depends on the value of the chemical po- 
tential. This is qualitatively completely different from HgTe, 
where the conductance displays a monotonic dependence on 
L. The oscillatory features become more smeared at elevated 
temperatures, as expected. 

Finally, we note that after this paper was submitted for pub- 
lication, we learned about the very recent work of H.-Z. Lu et 
a/J22l and C.-X. Liu et al!^\ in which similar finite size effects 
have been predicted. 



IV. SUMMARY 

In summary, we have investigated finite size effects on the 
surface states in the strong topological insulator Bi2Se3, com- 
paring the results also with those recently reported for Hg 
We demonstrate that the surface states respond differently to 
finite size effects in these materials, both quantitatively and 
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qualitatively. First of all, while the edge states become mea- 
surably gapped around L ~ 200 nm in HgTe^l the surface- 
states in Bi 2 Se 3 display a considerable robustness towards de- 
creasing L and become measurably gapped around L ~ 10 
nm. In this way, the topological surface state remains intact 
for a wider range of widths L. Secondly, the gapping between 
the surface states features a qualitatively distinct dependence 
on L in Bi2Se3 compared to HgTe, namely an oscillatory de- 
cay with L which stems from the material parameters that give 
an eigenvalue pair A Q that are complex conjugates. Both of 
these effects can be experimentally detected in a two-terminal 



geometry by varying the width and the chemical potential of 
the junction. 
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